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Abstract 


The aim of this paper is to study the dependence of reflection, absorption, and transmission of light by 
optically thick, weakly absorbing disperse media on the shape of constituent particles. Single scattering by 
nonspherical particles is calculated in the framework of the Rayleigh and geometrical optics approximations. 
The ratio of the average volume to the average surface of particles and their volumetric concentration are 
kept constant. Thus, the extinction coefficient and the optical thickness of different media under considera- 
tion are the same. Calculations are performed for ice particles at wavelengths of 0.645, 1.24, 1.64, and 2.13 um. 
The respective refractive indices range from 1.3082-i1.325 x 1078 (weak absorption) to 1.2674-i5.65 x 1074 
(moderate absorption). The radiative characteristics are calculated in the framework of approximate 
analytical formulae valid for optically thick, weakly absorbing media. © 1999 Elsevier Science Ltd. All 
rights reserved. 


1. Introduction 


Studies of the influence of the shape of aerosol and cloud particles on their optical properties are 
important for many applications in climate modelling [1], remote sensing [2,3], optical particle 
sizing [4], and spectroscopy of disperse media [5]. It follows from general considerations that this 
influence can often be neglected for particles that are much smaller than the wavelength of the 
incident radiation (a < 4/|m|, where a is the particle size, 2 is the wavelength, and m = n — iy is the 
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particle refractive index). However, this is not the case for large (a > 2) scatterers. For instance, 
Takano and Liou [6] demonstrated that neither area-equivalent nor volume-equivalent spheres 
can adequately approximate large hexagonal cylinders in terms of their single scattering and 
absorption properties. The same is true for circular cylinders and spheroids [3,7]. 

Differences in radiative characteristics of multiple-scattering media consisting of spherical and 
nonspherical particles are less pronounced than for single-scattering characteristics. This is due to 
the randomization of photon directions in the course of many scattering events. The investigation 
of the influence of crystal shapes on the radiative fluxes at visible wavelengths was performed by 
many authors [8-12]. Numerical solutions of the radiative transfer equation were used in these 
studies. It was found that the differences in radiative fluxes of media consisting of spherical and 
nonspherical particles decrease with increasing optical thickness t. 

The aim of this paper is to investigate the influence of particle nonsphericity on radiative fluxes in 
plane-parallel layers with large optical thicknesses (t > 7). This allows us to use simple asymptotic 
equations for the radiative characteristics [13,14]. We consider Rayleigh (a < 2) and large (a > A) 
particles. The following shapes of randomly oriented scatterers are studied in the framework of the 
ray optics approximation [15,16]: spheres (s), prolate spheroids (ps), oblate spheroids (os), circular 
columns (cc), circular plates (cp), hexagonal columns (hc), and hexagonal plates (hp). The ratio of 
the volume to the surface area of the particles is kept constant. Calculations for ice particles are 
performed at wavelengths 2 = 0.645, 1.24, 1.64, and 2.13 um. 


2. Single-scattering characteristics 
2.1. Rayleigh scatterers 


For Rayleigh scatterers (a < 2/|m|), the absorption cross section C,,, is proportional to the 
particle volume V (C,,, = «V) and the scattering cross section Cea is proportional to the particle 
squared volume V? (C,., = BV”). Furthermore, the phase function p(0) = 0.75(1 + cos? 0) does not 
depend on the particle shape and size. In the case of soft particles (m — 1), the dependence of the 
constants « and f on the particle shape is negligible. In general, x and p depend on the shape of the 
particles and on their refractive index m [17]. For example, for randomly oriented ellipsoids we 
have [17-19]: 


a =fi(n)y, (1) 
B = fr(n)k*, (2) 
where 
ne 1 
fe)=3 > apa ean 3) 
Salm) 1 3 w+ +DL + Ww? (4) 


~ 180, [(l + ob)? + Ll??? ’ 
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and v=n° — 1 — x°, w= 2ny, k = 27/2, y = 4ny/d. L; are called geometrical factors [17]. The 
absorption coefficient y is small for many substances in the visible region of the electromagnetic 
spectrum and Eqs. (3) and (4) simplify in this particular case as well: 


wou vest: 


AM=5 187 


fan) = y, (5) 


where 
3 
y= X [+ 0- DEI. (6) 


The relative errors A, and A, of the absorption and scattering cross sections due to neglecting the 
nonspherical shape of particles, 


ffi BP 
A, = ; A, = 7 
coincide in this particular case. They are equal to (see Eq. (6)) 
y” — y5 
A = ps (8) 
ys 


where indices n and s denote quantities for nonspherical and spherical shapes, respectively. For 
instance, for prolate spheroids with semi-axes r, b, and / (r = b < 1), it follows [19] that 


ie 1l+e 1 1-L, 1— Li 

L =f14 ye Ins, >| L, = 7” Lyss (9) 
where e = ./1 — &* and é = r/l. For small values of ¢ we have: 

1+e E g 

mofe 5E o) (10) 
and 

| 2e? 
Li ==- —. 11 
HgS (11) 


The geometrical factors for spheres (e = 0) are L; = L, = L} = 4 . Thus, L, for spheroids is 
always smaller than that for spheres. For a given value of ¢, the geometrical factors for oblate 
spheroids (and, thus, the light scattering characteristics) are closer to the respective values for 
spherical particles than those for prolate spheroids [19]. 

For soft particles (m > 1) at y/n < 1, it follows that fi(n) > 1, f2(n) > (2/37)|m — 1|? regardless of 
their shape. In this particular case only the particle volume is important. This is approximately true 
even for comparatively large refractive indices such as those of ice (n = 1.31) and dust (n = 1.52) at 
visible wavelengths (see Fig. 1). Thus, single light scattering characteristics of small nonspherical 
(a < A/|m|) particles can be approximately represented by those of volume-equivalent spheres. It 
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Fig. 1. Dependence of the relative error A [see Eqs. (8), (16), and (9)] on the aspect ratio č for Rayleigh ice (n = 1.31) and 
dust (n = 1.53) prolate spheroids in random orientation. 


appears that the same is true [1] even for larger particles (a < 54), but only with respect to the 
integral light scattering characteristics, namely the single scattering albedo wọ = 1 — (Cavs/Cext) 
(Gaps and Cex are the extinction and absorption coefficients), the extinction coefficient Cex and 
the asymmetry parameter g = 4{§p(0)sin 0cos0d0. According to the similarity principle [13], 
the radiative fluxes and albedos of scattering media depend mostly on the parameters 


s= J (1 — w,)/3(1 — g) and (1 — g). Thus, volume-equivalent spheres can be used to model the 
radiative fluxes of scattering layers with small to moderately sized particles. Of course, such 
a substitution is not possible if one needs to estimate polarization characteristics or bidirectional 
reflectances of a layer consisting of nonspherical particles [8-12,20]. 


2.2. Large convex particles 


The restrictions on the particle size and refractive index that will be considered now are as 
follows: x > 1, 2x|m — 1| > 1, where x = 2na,,/2, am is the minimal dimension of a convex particle. 
It is assumed that particles are monodisperse and randomly oriented. The extinction coefficient 
Oext for such large convex particles in random orientation can be found from the following 
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Table 1 
Volume V, surface area £, and effective radius a.t = 3V/X for particles of different shapes 


Type of a particle V 2 def 
Sphere (a - radius) 4na?/3 4ra? a 
Prolate spheroid 2] 2| 1, aresin,/1 — ¢ a 
(47/3)r 4nr] -+ 7 5 z 
(b =r < l are semi-axes, 2 2/1- ¢? È + (arcsin / 1 — 4/1 — ¢°) 


e=/1—@, €=7/l) 


In(l + e) — e’) r 
Oblate spheroid (4n/3)r71 in| + n=] ; ; 
ae IMs A=) 


2ée 


(b =r > l are semi-axes, 
&=/—1, €=r/l) 


Circular cylinder (/ is the nr’l 2ri(1 + né) 3nr/2(1 + né) 
length, r is the radius, 


€= r/) 


Hexagonal cylinder (I is J 3r?l/2 6rl(1 + (€/2,/3)) V3r 
the length, r is the side of 4(1 + (€/2,/3)) 


the cross section, € = r/l) 


equation [7]: 


fe (12) 


where aet = 3V/È is the effective radius, Vis the volume of a particle, X is the surface area, and c is 
the volumetric concentration. For spheres with radius a, the value of a,; is equal to a. Thus, it 
follows that in the case of large nonspherical particles the appropriate effective radius (at least for 
the calculation of the extinction coefficient and optical thickness t = Gextl, where /is the geometrical 
depth) could be the radius defined via the ratio of the average volume and surface area of particles. 
Extinction coefficients of spherical and large nonspherical convex particles in random orientation 
with the same values of a,, and c coincide. 

To check the applicability of a common effective radius to other light scattering character- 
istics, we have calculated the single scattering albedo wo, the asymmetry parameter g, and the 
similarity parameter s = J (1 — @o)/3(1 — g) for particles of different shapes by means of the 
geometrical optics approximation [15,16]. The ratio V/X was kept constant (aet = 100 um). 
The types of particles together with the correspondent definitions of V, X, and a.t =3V/X 
are presented in Table 1. The corresponding geometrical dimensions are shown in Table 2. 
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Table 2 

Parameters of particles used in calculations (b = r and l are semi-axes of spheroids, l is the length of 
a hexagonal or a circular cylinder, r is the side of the hexagonal cross section or the radius of 
a circular cylinder) 


Type of a particle Sizes The shape parameter 
Sphere (S) a = 100 pm 1 
Prolate spheroid (PS) r = 85.46 um, / = 170.92 um 2 
Oblate spheroid (OS) r = 138.02 um, / = 69.01 um 2 
Hexagonal column (HC) r = 93.65 um, | = 374.59 um 4 
Hexagonal plate (HP) r = 143.65 um, l = 143.65 um 1 
Circular column (CC) r = 83.33 um, / = 333.33 um 4 
Circular plate (CP) r = 133.33 um, / = 133.33 um 1 
Table 3 


Wavelengths 4, refractive indices of ice m = n — iy, and attenu- 
ation parameters p used in calculations 


A, pm n x P 

0.645 1.3082 1.325e — 8 5.1629e — 5 
1.24 1.2972 1.220e — 5 2.4727e — 2 
1.64 1.2881 2.670e — 4 4.0917e — 1 
2.13 1.2674 5.650e — 4 6.6667e — 1 


The shape parameter in Table 2 is defined as the ratio of the largest to the smallest particle 
dimensions. 

The wavelengths, refractive indices m = n — iy, and attenuation parameters p = 2a,¢y with 
y = 4ny/2, de = 100 um are presented in Table 3. The wavelengths selected are almost free of 
gaseous absorption and correspond to channels 1, 5, 6, and 7 of the Moderate Resolution Imaging 
Spectroradiometer (MODIS) on board of the NASA’s Earth Observing System satellites. 

One can see that the attenuation parameter p is smaller than 1 in our calculations. Large 
absorption is not of interest for this study because in this case wọ and g for spheres and randomly 
oriented nonspherical convex particles coincide [17]. The geometric optics results for the co- 
asymmetry parameter 1 — g, the co-single scattering albedo 1 — wo, and the similarity parameter 
s are presented in Table 4. 

First of all it follows that 1 — g and 1 — œọ differ considerably for spherical and nonspherical 
particles with the same effective radius a.r. Thus, in agreement with previous work on this subject 
[8-12], it is not possible to substitute nonspherical particles by V/X-equivalent spheres for 
calculations of the single scattering characteristics. This was demonstrated in our earlier papers 
[3,7] as well. For instance, for hexagonal plates at 2 = 0.645 um, 1 — g is about two times larger 
than that for spheres. The probability of photon absorption by spheres b = 1 — wo is about 1.5 
times smaller than that by hexagonal columns. 
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Table 4 


Dependence of values of 1 — g, 1 — wo, and s = ./(1 — wo)/3(1 — g) on the type of particles 
(see Table 2) for different wavelengths 2 and ast = 100 pm 


Shape 2, ym 1—g 1— wo s 

S 0.645 1.069 E — 1 2.158 E—5 8.203 E — 3 
PS 0.645 1.340 E — 1 2.497 E—5 7.888 E — 3 
OS 0.645 1.718 E-1 2.611 E—5 7.118 E—3 
HC 0.645 1.861 E-1 3.028 E—5 7.364 E — 3 
HP 0.645 2.004 E—1 2.921 E—5 6.969 E — 3 
CC 0.645 1.603 E — 1 2.646 E —5 7418 E—3 
CP 0.645 1.738 E — 1 2.635 E—5 7.108 E — 3 
S 1.24 1.016 E— 1 1.016 E—2 1.827 E-—1 
PS 1.24 1.261 E-1 1.164 E—2 1.754E-—1 
OS 1.24 1610 E-1 1.214E—2 1.585 E—1 
HC 1.24 1.743 E — 1 1.398 E— 2 1.635E— 1 
HP 1.24 1.902 E — 1 1.349 E- 2 1.537E— 1 
CC 1.24 1517 E — 1 1.237 E—2 1.649 E—1 
CP 1.24 1.660 E — 1 1.233 E—2 1.574E-1 
S 1.64 8.045 E — 2 1.398 E— 1 7.610 E—1 
PS 1.64 8.678 E — 2 1.507 E— 1 7.609 E—1 
OS 1.64 1.098 E — 1 1.559 E— 1 6.880 E — 1 
HC 1.64 1.166 E — 1 1.696 E— 1 6.963 E— 1 
HP 1.64 1.355 E — 1 1.671 E— 1 6.410 E—1 
CC 1.64 1.084 E — 1 1.584 E— 1 6.980 E— 1 
CP 1.64 1.246 E — 1 1.609 E — 1 6.560 E — 1 
S 2.13 6.584 E — 2 2.056 E— 1 1.024E+0 
PS 2.13 6.593 E —2 2.517 E-1 1.044E+0 
OS 2.13 8.105 E—2 2.215E-1 9.545 E—1 
HC 2.13 8.833 E—2 2.342 E-1 9401 E-1 
HP 2.13 1.066 E — 1 2.338 E— 1 8.549E — 1 
CC 2.13 8.424 E — 2 2.246 E— 1 9.427 E—1 
CP 2.13 9.828 E—2 2.289 E—1 8.810 E—1 


For all cases used in this study, 1 — g and 1 — œo are the smallest for spherical particles. 
Therefore, spheres are characterized by smaller absorption and less isotropic phase functions in 
comparison with V/X-equivalent convex nonspherical particles (at least for shape parameters 
n < 4; see Table 2). This is not the case for extremely elongated or flattened particles [7]. However, 
small to moderate deviations from sphericity (which are widely encountered) should produce this 
result in general. The largest differences in absorption are found between spherical and hexagonal 
particles, the smallest differences are between spheres and spheroids. Cylinders occupy the inter- 
mediate place. 

The values of 1 — g and 1 — wo vary considerably for different particle shapes at aet = const. 
However, the variation of the similarity parameter s = < (1 — œ0)/3(1 — g) is small (see Table 4). It 
remains the same within 10% for all shapes considered (at 2 = const.). The similarity parameter 
determines the spherical albedo of semi-infinite turbid media [13,14]. The larger values of s for 
spheres (see Table 4) imply smaller albedos. 
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3. Radiative characteristics of optically thick layers 


In order to study the influence of different particle shapes on the radiative characteristics of 
turbid layers, we make use of the single scattering results shown in Table 4. The radiative 
characteristics are calculated with approximate equations for weakly absorbing (s < 0.2), optically 
thick (t > 7) layers. These equations are derived in the framework of the radiative transfer theory 
[13,14,21]. They are presented in Table 5 for the plane albedo R(uo), the spherical albedo R, the 
transmittance T(t), the global transmittance T, and the global absorptance A = 1 — R — T. Here 
Ho is the cosine of the solar zenith angle. The accuracy of the formulae in Table 5 was discussed by 
Zege et al. [21]. The functions in Table 5 are defined as follows: 


Rito) == | dnf Roen wav, (13) 
R= i | pdn | po do) Ros Ho, Y) dy, (14) 
Tuo) == | ndu| Teno way, (15) 
T= z fs an| no do| Tos Ho, W) dy, (16) 


where u = cos & Y is the observation azimuth angle, 9 is the observation zenith angle, R(u, Uo, Y) is 
the reflection function of the layer, and T(u, Ho, Y) is the transmission function of the layer [13,14]. 
For nonabsorbing particles, the equations of Table 5 can be simplified: 


R(uo) = 1 — T (uo), T(Uo) = TK (Uo), (17) 


Table 5 
Radiative characteristics of weakly absorbing, optically thick layers [14,21] 


(y= 4s, s=4,/(L — @o)/3(1—g), 2= 31 —g)t, Ko(Mo) = 4(1 + 2H0)) 


Radiative characteristic Formula 
Absorptance A sinh(y ) + sinhQz2) 
sinh((1 + z)y) 

Plane albedo R(uo) exp( — yKo(Ho)) — TKo(Ho) exp( — (1 + z)y) 
Spherical albedo R — 

sinh((1 + z)y) 
Transmittance T (uo) T Ko(lo) 

sinh(y) 


Global transmittance T 5 ae 
sinh((1 + z)y) 
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and 
Z 1 
R= T = 1 
1+7 1+2 ve) 
where 
z= (1 — g). (19) 


According to Table 4, spherical particles have the smallest values of 1 — g (and z). Thus, the global 
transmittance of light by layers consisting of nonabsorbing nonspherical particles T, is lower [see 
Eq. (18)] than the global transmittance T, in the case of layers with spherical particles assuming the 
same t. The relative differences are 


T; = Ta _ 0.75t(9s > Gn) 


e aomi eM) 


where g, and g, denote the asymmetry parameters for spherical and nonspherical particles, 
respectively. Correspondingly, the reflection of light by nonabsorbing layers consisting of non- 
spherical particles R, is greater than the reflection R, by layers of spheres: 


R,— Rp | Jn — Js 


R +0751 — gla) (21) 


This qualitative behaviour also holds for weakly absorbing semi-infinite light scattering media. 
From Table 5 as t > œ, it follows that 


Rx(Uo) = exp( = VK o(Mo)), Ry = exp( = y), y= 4s, (22) 


where R,, and R,,(o) are the spherical and the plane albedos of a semi-infinite layer. As expected, 
the transmittance is equal to zero at t > œ. Again, with y = 4s being the largest for spheres (see 
Table 4), it follows that Rs,» < Ra,» where symbols n and s denote nonspherical and spherical 
scatterers, respectively. 

Increasing the particle nonsphericity produces the same effect on the radiative fluxes as 
decreasing their size. Smaller particles imply smaller absorption and larger reflection. For the 
solution of the inverse problem in the case of disperse media consisting of weakly absorbing 
nonspherical particles (e.g., snow [22-24]) it follows that the approximation of spherical particles 
will underestimate the particle size. This underestimation is often related to the presence of 
pollutants in snow (e.g., soot, dust, heavy metals, etc.). It is interesting to note that the parameter 
y is proportional to both the degree of particle nonsphericity and the concentration of pollutants. 
Indeed, the value of y is determined by 


1 — o 
3(1 = g)’ 


where 1 — Wo = Gap;/Cext. The extinction coefficient is determined by Eq. (12). The absorption 
coefficient of a large, weakly absorbing particle can be approximated by o,,, = fyc, where y = 4ny/A 


y=4 (23) 
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and fis the shape parameter [7]. We then have 


y = N Yaet, (24) 


where the parameter 


eee 
Le er (25) 


depends on the particle shape. Thus, the accuracy of the values of aep retrieved from the value of 
y = In(R;,') depends, according to Eq. (24), on the assumed value of N, /y: 


2 
det = ra (26) 
Therefore, low values of a.f as compared to in situ measurements may be explained by the fact that 
the true absorption coefficient and/or the true value of N are larger than those used in the retrieval 
procedure. 
It follows from Table 4 that N is always larger for nonspherical particles than for spheres. Note, 
that this important parameter can be found from Eq. (26) if the effective radius is available from 
independent in situ measurements: 


y 
~V Yaet 


Let us consider the absorptance A = 1 — R — T now. It follows for small values of y (see Table 5) 
that 


N= 


(27) 


A=, (28) 


where T, = Gapsl is the optical thickness due to absorption of light by particles. The value of 
Cavs = (1 — Wo) ext is smaller for spheres (see Table 4). Thus, the absorptance of light by layers with 
nonspherical particles is higher when y > 0. 

Figs. 2 and 3 show the results for R, T, and A at 2 = 0.645 and 1.24 um based on the formulas 
presented in Table 5 and the data from Table 4. As expected [12], the differences in T(t) and R(t) 
for different shapes generally decrease with increasing optical thickness, being smaller for smaller 
absorption. This can be explained by the fact that photons in semi-infinite media are highly 
randomized and the shape of particles is not so important as in the case of thinner layers. 

The absorptance A is most sensitive to the particle shape. Layers with spherical particles absorb 
less radiation than layers with nonspherical particles at small y. This is related to the fact that A is 
proportional to the absorption coefficient cps at small values of y [see Eq. (28)]. However, for 
t> ©, we have [see Eq. (22)] 


A=1-—exp(-— y). (29) 


It follows that regardless of the fact that the absorption coefficient of spherical scatterers is 
smaller, the absorptance of the radiation in a semi-infinite layer with spherical particles 
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Fig. 2. Dependence of the spherical albedo (a), global transmittance (b), and absorptance (c) for layers with different types 
of particles (s — spheres, ps — prolate spheroids, os — oblate spheroids, hc — hexagonal cylinders, hp — hexagonal 
plates, cc — circular cylinders, cp — cicular plates) on the optical thickness at a wavelength of 1 = 0.645 um. 
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Fig. 3. As in Fig. 2, but for 2 = 1.24 um. 
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is higher than in case of layers with nonspherical scatterers. This is because y = 4s for spheres is the 
largest among all particles under consideration (see Table 4). Thus, semi-infinite layers with 
spherical particles absorb more radiation than semi-infinite media with nonspherical particles 
having the same effective radius. This result is due to the larger asymmetry parameter for spheres, 
which allows photons to penetrate into deeper layers and be trapped there. 

Opposite signs of the albedo difference = = A, — A, (A, and A, are the absorptances of layers 
with spherical and nonspherical particles, respectively), at small and large thicknesses imply that 
there is a certain value of t = tọ where 5 = 0. Note that to depends on the particle shape (see Figs. 
2c and 3c). 

It follows from Figs. 2 and 3 that weakly absorbing thick layers with spherical particles appear 
darker (smaller values of R) than layers with nonspherical particles with the same volume/surface 
ratio and concentration. They are also more transparent (larger values of 7). This agrees with 
results for nonabsorbing layers [see Eqs. (20) and (21)]. Results for elongated and flattened particles 
of the same type (cylinders, hexagonal particles) are very close. In all cases the radiative fluxes for 
prolate spheroids are closest to those for spheres. On the other hand, differences in R, T and A are 
most pronounced for spheres and hexagonal columns. 


4. Conclusion 


It is a well-known fact that bidirectional reflection and transmission functions of turbid layers 
consisting of nonspherical particles strongly depend on their shape [8-12,25]. The same is true for 
their polarization characteristics. 

However, it is sometimes believed that the difference between integral radiative character- 
istics for particles of different shapes is averaged out by multiple scattering. This is indeed the 
case for small (compared to the wavelength of the incident radiation) particles (see Fig. 1). 
In this case, nonspherical particles can be substituted by volume-equivalent spheres. However, 
for large particles the difference between different shapes of scatterers must be taken into 
account. 

We investigated the special case of weakly absorbing, optically thick (t > 7) layers with particles 
of different shapes. It was found that semi-infinite, weakly absorbing turbid media consisting of 
spherical particles are characterized by smaller spherical albedos R and larger absorptances 
A compared to layers with nonspherical particles with shape parameters smaller than 4 and 
der = const. It is interesting to note that these differences are smaller for randomly oriented infinite 
circular cylinders [26]. 

The absorptance of light at small values of t, (y > 0) is smaller in the case of media consis- 
ting of spherical particles (see Figs. 2c and 3c). The transmittance of light through finite layers 
with spherical particles is higher than in the case of layers with nonspherical particles. 
Taking into account different incidence angles (see Table 5) does not change our con- 
clusions. 

Differences between phase functions, single scattering albedos, and asymmetry parameters of 
spheres and randomly oriented convex nonspherical particles gradually disappear with increasing 
absorption [17,27]. In this limit the model of equivalent spheres is indeed appropriate. This factor 
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simplifies all calculations at large absorption and can be used, e.g., in the spectroscopy of turbid 
media with large, strongly absorbing nonspherical particles [5]. 
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